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We propose a strong-coupling analysis of a polarizable planar interface, in the spirit of a recently introduced 
Wigner-Crystal formulation. The system is made up of two moieties: a semi-infinite medium (z < 0) with 
permittivity e' while the other half space in z > is occupied by a solution with permittivity e, and mobile 
counter-ions (no added electrolyte). The interface at 2 = bears a uniform surface charge. The counter-ion 
density profile is worked out explicitly for both repulsive and attractive dielectric image cases. 
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Introduction. Strong coulombic couplings can conveniently be realized in colloidal systems by increasing 
the valence of micro-ions HI 121 . The counter-intuitive effects such as like-charge attraction or overcharging 
that ensue |[T]|2][3], have been mostly studied under the simplifying assumption that the dielectric permittivity of 
colloidal particles coincides with that of the surrounding solvent -assumed to be a structure-less medium, charac- 
terized by its sole static dielectric response. However, realistic systems, of biological interest in particular, often 
exhibit highly inhomogeneous dielectric structure, with polarizable colloids having a static dielectric constant 
{e' < 10), much smaller than that of the solvent (e ~ 80 for water). It is therefore desirable to include colloidal 
polarizability into existing strong-coupling treatments (see below for a precise definition of the relevant coupling 
parameter). This is the purpose of the present contribution. 

To our knowledge, the dielectric inhomogeneities were studied in the leading strong-coupling order only for 
the geometry of two parallel plates ID where counter-ions are confined in a slab, and similarly in Ref. ID, for a 
spherical macro-ion in a concentric spherical Wigner-Seitz cell. It was realized that while the image charge effects 
are generally small in the weak-coupling regime, they become relevant in the strong-coupling limit. Our goal here 
is to obtain exact results for a single polarizable planar colloid (hatched region in Fig. [T]), in the strong-coupling 
limit (see Fig. [T]for a definition of the coupling parameter S). To this end, the Wigner-Crystal formulation of 
Strong-Coupling (SC) theory, recently introduced ||6l, is modified to account for dielectric image charges. We 
emphasize that this approach has been shown to be in excellent agreement with Monte Carlo simulations, in 
contrast to previous fugacity SC expansions lO. The idea is to consider the first relevant large S excitations 
around the ground state structure, that is a Wigner Crystal realized when, strictly speaking, S — > 00 @. 
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Fig. 1 The one-plate geometry considered. Counter-ions of charge — ge, e being 
the electron charge, are immersed in a solvent of dielectric constant e in the region 
z > while the wall, a large colloid with surface charge ere occupying the z < 
half-space, has the dielectric constant e' . From these quantities and the inverse tem- 
perature /?, it is customary to define the coupling parameter H = 2-Kq'^(.%a, where 
£s = /3e^/e is the so-called Bjerrum length. When H <C 1, mean-field treatments 
are legitimate, while the strong-coupling regime of interest here corresponds to the 
limit of large H. 
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We define the dielectric jump between the two moieties as A = (e — e')/(e + e'). The image of a charge e at 
position r = {x,y, z > 0) has the charge of strength e* = eA and is locaHzed at position r* = (x, y, —z < 0). 
Counter-ions interact via the Coulomb interaction potential u(r, r') = uo{r, r') + Mim(r, r'), where wo(r, r') and 
Wim(r, r') are the direct and image Coulomb potentials, respectively, defined by 



wo(r,r') = 
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The total interaction energy of the counter-ions at positions {rj}^^^ in the half-space z > reads 
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Here, the two one-body terms in the square brackets originate in the interaction of the particle with the surface 
charge plus its image (therefore the factor 1 + A) and in the interaction of the particle with its own charge image. 
We shall study separately the cases of attractive (A < 0) and repulsive (A > 0) image charges. In addition, 
we shall consider dimensionless quantities: z ~ z/ i^l where /i = 1/ {inqlBa) is the Gouy-Chapman length, and 
likewise, counter-ion densities will be made unit-less by the transformation p = pj {^-kIbo^). 

Attractive images. When A < 0, a small hard core of radius h is required to prevent the thermodynamic 
collapse of charges onto their images at the surface of the wall. Irrespective of this steric regularization, the 
ground state is the same as when A = 0. It corresponds to the usual hexagonal Wigner layer of counter-ions 
with the lattice spacing a given by the electroneutrality condition q/a ^ \fZa} jl, localized at the hard-core 
distance h from the wall |7 !|. We shift particles from their ground-state Wigner positions by small displacements 
Zj . Expanding the two-body potential in (|2l) up to harmonic particle shifts, forgetting irrelevant constant terms, 
and keeping in mind that z > 6 everywhere, the energy change from the ground state is given by 
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Here, all particles are indexed by an integer, and Rjk denotes the distance between ions j and k in the ground 
state structure. In the remainder, we will refer to the one body term on the r.h.s. of the first equation in Q, 
corresponding to the square brackets, as —f3SEo. In other words, —(3SE = —(35Eo + Sz- The associated one- 
body Boltzmann factor is dominant at large S, and reads cxp[— (1 + A)z + S|A|/(4z)]. Up to a normalization 
factor, it therefore provides the leading order term in the counter-ion density expansion, that we denote po{z). As 
expected with attractive images, this contribution stems from the region of z close to b. 

In order to compute the first sub-leading correction to the density profile, we add to the one body potential 
SEq a generating potential u(r), to be set to at the end of the calculation. We next resort to the cumulant tech- 
nique sketched in [6|, where, having introduced iu{r) = exp(— /?u(r)) the Boltzmann weight of the generating 
potential, the ionic density can be written as 



p(z) = poiz) + 
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In this relation, the dots are for higher order corrections, i.e. terms of larger power in and the statistical 
average (...)o refers to the system of non-interacting particles with Boltzmann weight exp(— /3(5i?o). Within this 
formalism, in the first correction order we have 
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The density profile then takes the form 
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where 
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This density profile satisfies the electroneutrality condition J p dz ^ a / q, which translates into / pdz ~ 1. For 
the homogeneous dielectric case A = and hard core 6 = 0, we recover the result obtained in |6l- For strictly 
negative values of A, the effect of attractive charge images is fundamental in the SC limit. 

Repulsive images. When A > (repulsive image charges), the counter-ions are on the one hand attracted to 
the wall by the surface charge tre and on the other hand repelled from the wall by their images. It is natural to 
assume that the ground state of the system is formed by the standard hexagonal Wigner crystal of counter-ions 
with lattice spacing a given as previously by the electroneutrality condition qja = \/ia? /2, localized at some 
distance / from the wall. This distance is determined by balancing the above attractive and repulsive forces. 

The energy per particle for the Wigner crystal at distance z from the waU is given by 
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On the rhs, the first term describes the interaction of the particle with the surface charge, the second (z-inde- 
pendent) one the direct interaction with all other particles, and the last two terms the interaction of the particle 
with its own self-image and with images of all other particles. The distance of the Wigner crystal from the wall 
I is determined by the stationarity condition dz\fiEQ{z) / N]\z=i = 0. Introducing the variable t = 2l/a, this 
requirement can be written as 
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where use was made of the hexagonal structure. The lattice sum appearing here can be represented in a more 



convenient way for numerical calculation. If j + A: is an even number, we introduce integers n 



and m = (j — k)/2 which "diagonalize" the expression + jk 
number, we introduce integers ^ {j + k + l)/2 and m ^ {j — k - 
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3(n - 1/2)2 + _ 1/2)2. Finally, using the formula T{s)[f{n, m)]-' = duM'-ie-^^^"'"), we find that 
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where 6*2(0, q) = Y.n g^""^/^^ and 6*3(0, q) = Y.n <?" ^re the Jacobi theta functions ID, the sums running from 
—00 to 00. The integral over theta functions can be evaluated with a high precision and very quickly, e.g. by 
using the symbolic language Mathematica. The numerical solution of Eq. ( fTOl ) for the extreme case A = 1 is 
t ~ 0.295. In the limit of small A, we have t ^ 3^/'^-\/A/(47r). For a non-negligible dielectric jump, with A of 
order unity, the distance of the Wigner crystal from the wall Z = at/2 is of the order of the Wigner lattice spacing 
a. 

We still need to examine the ground-state stability of the above Wigner crystal along all spatial directions. The 
stability in the {x, y) plane is related to the already made normal mode analysis in Ref. ||9]. To study the stability 
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of the Wigner crystal along the z direction, we shift one particle (say j = 1) from its lattice position by a small 
— z — I. The corresponding change in the energy 8E is given by 
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Expanding this expression in 6z, the linear term vanishes due to the stationarity condition (flOl l. The quadratic 
term reads as (36E = ((Jz)^/^^, where ^ is given by 
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The criterion for the thermodynamic stability of the Wigner structure is that the sign of the contributions in curly 
brackets be positive, and it is so in the whole interval A e [0, 1]. For large but finite S, the leading particle 
density profile is represented by the Gaussian distribution 



p{2 
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Since from (fTsT i the standard deviation ^ scales like ^/a oc 1 /S^/^, in the limit S 00 the Gaussian distribution 
becomes the (5-function peaked at z = / and the particles are constrained to the ground-state Wigner lattice as it 
should be. 

Conclusion. In this work, we have obtained the ionic profiles of counter-ions strongly coupled to a planar 
charged interface, that defines the boundary of a large colloid which has a static dielectric constant differing 
from that of the solvent (see Fig. [T]!. At variance with the Poisson-Boltzmann picture where dielectric images 
do not affect the profiles, the strong-coupling (large S) limit is highly sensitive to the dielectric jump A = 
{e — s')/{e + e') 7^ 0. When image charges are attractive (A < 0), a situation which requires the introduction 
of an ionic hard-core to avoid instabilities, the density profile is expectedly enhanced in the immediate vicinity 
of the wall (see Eq. d?)). On the other hand, with repulsive images (A > 0), the ground state of the system is 
made up of an hexagonal Wigner crystal at finite distance from the wall. Finite but large values of S turn the 
ground state configuration of counter-ions into a Gaussian form with standard deviation proportional to aj'E}!^, 
where a is the lattice spacing in the ground state Wigner crystal. It would be desirable to test quantitatively our 
predictions against numerical simulations, since those reported in a recent analysis f5l, in addition to pertaining 
to a curved geometry, also make use of a ion/wall hard core that prevents the study of the Gaussian behaviour 
embodied in Eq. ( fT4] i. The two plates configuration is also an interesting venue for future work. 
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